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Summary. We present an approximate, analytical formulae for the computation 
of hadronically-produced gamma-ray luminosity in the innermost regions of stellar 
winds. We put this into the context of other recent research on regions of star 
formation as gamma-ray sources. 

1 Introduction 

It has been proposed that several 7-ray sources are related with early-type 
stars and their neighborhoods, since collective effects of strong stellar winds 
as well as of supernova explosions are expected to result in particle acceler- 
ation up to energies in excess of a TeV (see, e.g., Montmerle 1979; Casse & 
Paul 1980; Bykov 2001, Romero & Torres 2003, Torres et al. 2003). Recently, 
several TeV unidentified source were detected spatially close to regions of star 
formation, most especially the first such source, detected in the Cygnus region 
(Aharonian et al. 2005), where a nearby EGRET source (3EG J2033-I-4118) 
has a likely stellar origin (White & Chen 1992; Chen et al. 1996; Romero 
et al. 1999; Benaglia et al. 2001). The TeV discoveries and the increase in 
sensitivity at hand in the GeV and TeV regimes has rekindled the interest 
in regions of star formation as gamma-ray sources. But detecting new popu- 
lations of gamma-ray sources is not a trivial task, particularly for surveying 
telescopes. For instance, a large fraction of the detections to be made by the 
Gamma-ray Large Area Space Telescope (GLAST) will initially be unidenti- 
fied, and, due to their sheer number, traditional methodological approaches to 
identify individuals and/or populations of 7-ray sources will encounter proce- 
dural limitations. These limitations will hamper our ability to classify source 
populations lying in the anticipated dataset with the required degree of con- 
fidence, particularly those for which no member has yet been convincingly 
detected in the predecessor experiment EGRET. A new paradigm for achiev- 
ing the classification of 7-ray source populations based on the implementation 
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of an a priori protocol to search for theoretically-motivated candidate sources 
is probably needed in order to protect the discovery potential of the sam- 
ple (Torres & Reimer 2005). Key to the any new procedure is a statistical 
assessment by which the discovery of a new population can be claimed un- 
der a similar quantitative technique. For the previously cited protocol, this 
is achieved by introducing that the hypothesis to test in searching for new 
population of sources is the null hypothesis against a reduced noise level. It 
is useful to note that the LAT in GLAST will be in a privileged position to 
actually identify new population of sources. If LAT would have an additional 
order of magnitude better sensitivity, with no very significant improvement 
in angular resolution, a situation similar to the GRB case would appear, i.e., 
a flat distribution of unidentified sources with a few privileged individuals 
only which are extensively traced in multifrequency studies. Essentially, we 
would find a 7-ray source coinciding with the position of every member of 
any population under consideration. And thus, we would lack the capability 
to achieve discoveries by correlation analysis. This is, perhaps, already indi- 
cating that a next generation high energy 7-ray mission after GLAST-LAT 
should not be exclusively sensitivity-driven if no significant improvement in 
angular resolution can be achieved. 



2 Hadronic processes in single and collective winds 

Regions that currently are or have been subject to a strong process of star 

formation are good candidates to be 7-ray and neutrino emitters and may 

even be sites where ultra high energy cosmic rays are produced (see Torres 

& Anchordoqui 2005 for a review). Outside the Galaxy, the more powerful 

sites of star formation are found within very active galaxies such as starbursts 

and Luminous or Ultra-Luminous Infrared Galaxies, and they have also been 

recently analyzed as high energy emitters (e.g. Torres et al. 2003, Torres 2004, 

Torres et al. 2004ab, CilUs et al. 2005, Domingo-Santamaria & Torres 2005). 

Here, we will not cover these results concerning collective effects in extragalac- 

tic scenarios, referring the reader to the mentioned publications for details. 

Rather, we will make a short note on the computation of hadronic gamma-rays 

in single and collective stellar winds. 

Typically, a stellar wind consists of an inner region in free expansion, 

followed by a larger, low dense region of hot shocked wind (bubble), which is 

finally surrounded by a thin layer of swept-up interstellar medium (ISM) (see, 

e.g., Lamers & Cassinelli 1999, Ch. 12). We assume that 7-ray production is 

mainly hadronic and that it occurs at the end of the free expansion phase of 

the wind evolution, when the swept-up matter is comparable with the mass in 

( ■ , \^/^ 
the wind driven wave. Its size is then i?wind = ( 3M*/47rmpnoVoo 1 , where 

Af* is the star mass loss rate, mp is the mass of the proton, no is the ISM 
number density and 14o is the wind terminal velocity. For a star of radius 
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Ri, the particle density and particle velocity radial profiles of the wind at 

this inner zone can be modeled by n(r) = (1 — Ro/r)~^ M^,/ A-KiripVaor'^ , and 
V(r) = 1^00 (1 — Ra/rY ; respectively. The parameter /3 is ~ 0.8 — 1 for massive 
stars, -Ro is given by i?o = — {Vq/VooY^^), and the wind velocity close to 
the star, Vq, is typically ~ 10~^yoo- As a first exercise, we can compute the 
resulting 7-ray luminosity of a single star with the target matter being that 
of the innermost region, and between energies Ex and E^: 

L^= I I n{r)q^{E^)E^{4-iTr'^)drdE^, (1) 

JR^ JEi 

with the differential 7-ray emissivity from Tr^-decays, q^{E^), approximated 

by 

q^{E^) = A7ra„{Ep){2Z^;'2^,/a) Jp{E^)r,AO{Ep - £™) (2) 

at the energies of interest. The parameter rjA takes into account the contribu- 
tion from different nuclei in the wind (for a standard composition t/a ~ 1-5, 
Dermer 1986). Jp{E^) is the proton flux distribution evaluated at E = E^ 
(units of protons per unit time, solid angle, energy-band, and area), which is 
assumed as a typical power-law Jp{Ep) = {c/4,iT)KpEp~°' , with a being the 
spectral slope and Kp the spectrum normalization factor. The cross section 
app{Ep) for pp interactions at energy E.p « lOii^-y can be represented above 
Ep K. 10 GcV by app{Ep) « 30 x [0.95 + 0.06 log(£;p/GeV)] mb (e.g., Aharo- 

nian & Atoyan 1996). Z^"^^,, is the so-called spectrum- weighted moment of 
the inclusive cross-section. Its value for different spectral indices a is given, for 
instance, by Drury et al. (1994). Finally 0{Ep — E'!^™) is a Heaviside function 
that approximately takes into account the fact that only CRs with energies 
higher than S™'"(r ^ i?*) will diffuse into the wind (see below). For normal- 
ization purposes, we take the expression of the proton energy density, which 
is given by 

WCR = j Np{Ep)EpdEp = {134.4 ; 9.9 ; 0.8}^!:^ eV cm'^ = ^wcr,©, (3) 

where <; is the enhancement factor of the CR energy density with respect to 
the local value, wcr,©, we have considered energies between 1 GeV and 20 
TeV, and the choice for the prefactor, {134.4; 9.9; 0.8}eV cm~^ , refers to 
three different slopes, 1.9, 2.0, and 2.1, respectively. If the prefactor is referred 
as f{a) [with units of a;cR,0, typically given in eV cm^'^J, the normalization 
is Kp{a) ~ LOcR,Q<i / f {ct) , depending on the slope of the spectrum. We assume 
the Earth-like spectrum to be Jq{E) = 2.2£'Qgy ^cm-2s~^sr-iGeV~\ so that 
i-^CR.0 ^ l.SeVcm^''. An analytic expression for the 7-ray luminosity of a 
single star produced in the innermost region of its wind can thus be obtained 
if the radial dependence of the wind velocity is ignored. This, in fact, does not 
introduce significant cliangc;s on the results because a) the terminal velocity 
of the wind is reached in only a few star radii, and b) the integration extends 
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to radii up to i?wind ^ -R*) with the largest radii (where most of the mass is 
located) dominating the sum. The expression for the approximate luminosity 
then reduces to: 



9.6 X 10-38 cM*z(:2.o??Ai^p(2) 
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for a=2, and 



9.6 X 10-38 c z(^^o »?Ai^p(a) 

-i^7 ~ 77 (-flwind — R*) X 





for a 7^2, with the proton mass measured in kg, the terminal velocity and the 
speed of light measured in cm s^^, the mass loss rate measured in kg s^^, 
and the wind and stellar radius given in cm. Analogously, the expected 7-ray 
flux over an energy threshold i?th can be evaluated taking into account the 
distance D from the Earth to the stellar association: 



F^{E^ > Eth) = / / n(r) q^{E^) inr^ dr dE^. 



(5) 



Also here an analytic - although approximate- expression, can be obtained: 
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which is valid for all a 7^ 1. Here, D must be given in cm, whereas all other 
quantities must be given in the same units as those used in Eqs. (4) to respect 
consistency with the numerical factor. This formulae are useful to fix the scale: 
a typical single star would produce, by interaction between cosmic rays and 
target matter in the innermost region of its wind, a flux one to three orders of 
magnitude (depending on the kind of star) below the sensitivity achieved by 
current telescopes in the TeV regime. Of course, this estimation does not take 
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into account other processes, particularly those happening in binary systems 
where two strong stellar winds interact. 

In the case of a collection of winds, a hydrodynamical model taking into 
account the composition of different single stellar winds needs to be adopted. 
As we see below, difFcrcnccs with the single stellar wind case arc notable. 
One possibility is to consider that there are N stars in close proximity, uni- 
formly distributed within a radius Rc, with a number density n = SN/AttRI . 



the surrounding ISM and with each other, filling the intra-cluster volume 
with a hot, shocked, collective stellar wind. A stationary flow in which mass 
and energy provided by stellar winds escape through the outer boundary of 
the cluster is established. For an arbitrary distance R from the center of the 
cluster, mass and energy conservation imply 



where p and V are the mean density and velocity of the cluster wind flow 
at position R and h is its specific enthalpy (sum of internal energy plus the 
pressure times the volume). Solving these equations with the corresponding 
boundary conditions, one obtains the density and velocity of the cluster as 
a function of distance from its center. This is presented in detail in the pa- 
per by Domingo-Santamaria & Torres (2006). Again, it is to be noted that 
not all cosmic rays will be able to enter the collective wind. The difference 
between an inactive target, as that provided by matter in the ISM, and an 
active or expanding target, as that provided by matter in a single or a col- 
lective stellar wind, is indeed given by modulation effects. The cosmic ray 
penetration into the jet outflow depends on the parameter e ~ VR/D, where 
V is velocity of wind, and D is the diffusion coefficient, e measures the ratio 
between the diffusive and the convective timescale of the particles. Taking into 
account such effects, we have studicxl collective wind configurations produced 
by a number of massive stars, and obtained densities and expansion veloci- 
ties of the stellar wind gas that is target for hadronic interactions in several 
examples (sec details in Domingo- S ant aman'a & Torres 2006). We have com- 
puted secondary particle production, electrons and positrons from charged 
pion decay, electrons from knock-on interactions, and solve the appropriate 
diffusion-loss equation with ionization, synchrotron, brcmsstrahlung, inverse 
Compton and expansion losses to obtain expected 7-ray emission from these 
regions, including in an approximate way the effect of cosmic ray modulation 
(see Torres & Domingo-Santamari'a 2005 for details on the code). Examples 
where different stellar configurations can produce sources for GLAST and the 
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MAGIC/HESS/ VERITAS telescopes in non-uniform ways, i.e., with or with- 
out the corresponding counterparts were found. Cygnus OB 2 and Westerlund 
1 maybe two associations where this scenario could be at work 
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